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Abstract 

Recent investigations seem to favor a cosmological dynamics according to which the acceler- 
ated expansion of the Universe may have already peaked and is now slowing down again As 
a consequence, the cosmic acceleration may be a transient phenomenon. We investigate a toy 
model that reproduces such a background behavior as the result of a time-dependent coupling 
in the dark sector which implies a cancelation of the "bare" cosmological constant. With the 
help of a statistical analysis of Supernova Type la (SNIa) data we demonstrate that for a certain 
parameter combination a transient accelerating phase emerges as a pure interaction effect. 
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I. INTRODUCTION 



By now there is large direct and indirect evidence that the Universe entered a period 
of accelerated expansion. Direct evidence is provided by the luminosity-distance data of 
supernovae of type la (see, however, {3]), indirect evidence comes from the anisotropy 
spectrum of the cosmic microwave background radiation (CMBR) [4], from large-scale- 
structure data [5], from the integrated Sachs- Wolfe effect p, from baryonic acoustic 
oscillations {?| and from gravitational lensing According to the currently prevailing 
interpretation, based on Einstein's GR, our Universe is dynamically dominated by two so 
far unknown components, dark matter (DM) and dark energy (DE). The latter contributes 
roughly 75% to the total energy budget, the former about 20%. Only about 5% are in the 
form of conventional, baryonic matter. The preferred model is the ACDM model which 
also plays the role of a reference model for alternative approaches to the DE problem. 
While the ACDM model can describe most of the observations, there still remain puzzles 
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and references 



9|. The present situation in the field is characterized, e.g., in 
therein. 

Both DM and DE manifest themselves observationally only through their gravitational 
action. In most studies both these components are regarded as independent. However, 
there exists a line of investigation that explores the consequences of a coupling within the 
dark sector. Corresponding models show a richer dynamics than models with separately 
conserved energy-momentum tensors. It has been argued that ignoring a potentially ex- 
isting interaction between both components may give rise to misinterpretations of obser- 
vational data regarding the equation-of-state (EoS) parameter 13j. A coupling between 
DM and DE may also be relevant with respect to the coincidence problem (see, e.g. [jjj 
and references therein). Models with an interaction matter-dark energy were introduced 
by Wetterich 15) . Meanwhile there exists a still growing body of literature on the subject 



-see, e.g. 



161 ] and references therein. Limits for the admissible interaction strength have 



been obtained for various configurations 17| • On the other hand, the transition from de- 
celerated to accelerated expansion can be understood as a pure interaction phenomenon 
in the context of holographic DE models 181-120 1. 
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Recently it was argued that the latest SNIa could favor a scenario in which the cosmic 

n 

acceleration has past a maximum value and is now slowing down again [1[ . In such a case 
the accelerated expansion could be a transient phenomenon with a late-time dynamics, 
incompatible with that of the ACDM model. It is this possibility which we are going 
to study in the present paper. We mention that models of transient acceleration were 



already discussed in 2jj and more recently in 22]. 

It is the purpose of this paper to provide a toy model which describes a transient 
cosmological acceleration as the consequence of an interaction between dark matter and 
dark energy. Such a dynamics cannot be obtained if the interaction represents a small 
correction to, say, the ACDM model. For models of this type the long-time cosmological 
dynamics will always be determined by the cosmological term and result in accelerated 
expansion. To achieve transient accelerated expansion, a twofold role of the interaction 
is necessary. At first, it has to cancel the "bare" cosmological constant and at second 
it has to generate a phase of accelerated expansion by itself. Acceleration has to be an 
interaction phenomenon. We show that these requirements can be fulfilled by interaction 
terms that combine powers and exponentials of the cosmic scale factor. 

The paper is organized as follows. In Sec. [Ill we introduce the basic interacting fluid 
model and study the influence of the coupling, described by an initially arbitrary function 
of the scale factor, on the cosmic acceleration equation. Sec. II I II investigates two interac- 
tion types between DM and DE, each of them characterized by two parameters, in detail. 
Conditions for the existence of transient acceleration are found under the condition of a 
vanishing total effective cosmolo gica l constant. In Sec. II VI we use a statistical analysis of 



the gold sample of the SNIa data 23] to find the admissible parameter ranges for transient 
cosmic acceleration. The best-fit models for both interactions are shown to be compatible 
with the ACDM model. Our conclusions are summarized in Sec. [V] 

II. INTERACTING FLUID DYNAMICS 

We model the present cosmic substratum as a mixture of three components: dark 
energy (DE) with an energy density p x , dark matter (DM) with an energy density p m and 
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baryonic matter with an energy density p&. The Friedmann equation in a spatially flat 
universe is 

3H 2 = 8nG ( Px + p m + Pb ) , (1) 

where H = a /a is the Hubble parameter and a is the scale factor of the Robertson- 
Walker metric. The DE component is supposed to be characterized by an equation of 
state p x = wp x . We admit an interaction between DE and DM according to 

p x + 3H{1 + w)p x = -Q , 

p m + 3Hp m = Q , (2) 

where Q symbolizes the interaction, and a dot represents a derivative with respect to the 
cosmic time. The baryon component is separately conserved, 

p b + 3Hp b = pb = Pb a (— ) • (3) 

The split of the total energy density p into p = p m + p^ + p x is consistent with the energy 
conservation for p. We assume, without loss of generality, that the energy density of the 
dust fluid can be written as: 
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Pm Pmo 



(?) ><•> • w 



where p mo and a are constants, and f(a) is an arbitrary time-dependent function. It 
follows from P) and flU) that 

,— = n, — — 



Q = Pmj = Pm { — ) f ■ (5) 

Let us write f(a) as 

f(a) = l+g(a) , (6) 

such, that the deviation from the non-interacting case is encoded in the function g(a). 
With 

* ■ dg ■ (f\ 

f = 9 = Ta a (7) 



we obtain 

_ dg . (clq 
Q = Pm -j-a — 
da \ a 



For p m it follows that 



Pm 0-+9) 



(9) 



Denoting the value of p m at a = ao by p mo , both initial values are related by 



Pm = Pm (1 + 9o) 



(10) 



where go = g(ao). The quantity p mo is the value of p m at a = ao in the presence of the 
interaction, p mo is the value of p m at a = ao for vanishing interaction. The interaction 
re-normalizes the present value of p m . According to OH]), a positive value of Q, i.e., an 



energy transfer from dark energy to dark matter requires $ > in an expanding universe. 
For < the transfer is in the opposite direction, i.e., Q < 0. 

With (jBJ), the energy density of the x-component is determined by 

p x + 3H(l + w)p x -- 

For a constant EoS parameter w the solution of ( 1TTT) is 



Pmo j a I 
da V a 



Px = P. 



^Y il+W) -p m y o a-^ f da^a 3 " 



a J 



da 



(12) 



Integrating by parts yields 



px = (,Px + p mo go) 



a \ 3 ( l + w ) _, / a o\ 



a / 



Pm„ | — I 5 , + 3wp mo aoa 3(1+u,) / dagu 



3w-l 



(13) 



For a given interaction g(a), the Hubble rate in Eq. ([[]) is then determined by the sum of 
p b from <$\), p m from (j5J) and p b from f|T3|) . 
The acceleration equation 

d AirG 
a 

takes the form 

a 



a 
a 



\-c; f . ///..N--! 

i Pmo (1 + flO . ; 
3 V a / 



[Pm + Pb+ (1 + 3w)p a 



^Vp,^)%d + 3.) 



(14) 



(Pxo + Pmo^o) ( —J 



-p mo (^) 2 + 3«;p mo a 3 a- 3 ( 1+ ^ jf da^a 3 ^ 1 



(15) 
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It can also be written as 



a 
a 



— — Hn { fir 



a 



n 



bo 



«0 



+3w Ct % 



a 



+ (1 + 3w) 
1 + 3w) a~ 3w 



fix + ^m 90 

dag a 3 ™' 1 



a Q \ 3(i+io) 



where we have introduced the parameters Cl r 



8nGpm q 
3H$ ' iim o 



8iTGpm 



(16) 



87rGp3; 



and fibo = 8 ■ F° r an Y specific model g"(a), formula (jTB"]) describes the dynamics of 2 
through the cosmic history, as long as our three-component model makes sense. 
The present value of the deceleration parameter is 

1 



-A 



n 



mo 



[i + 9o ) + n bo + (i + 3w)n 



.Co 



{fi mo + tt bQ + (1 + 3w) flj-J 



(17) 



where Vt mo = Q mo (1 + g ). In a spatially flat universe with Vt mo + f2b = 1 — Q Xo , Eq. ([T 
reduces to 

^5 |o=-^{l + 3 W fi, } . (18) 
Notice that there is no direct influence of the interaction on the present value of the 



deceleration parameter 



24]. 



III. TRANSIENT ACCELERATING PHASES IN THE UNIVERSE 



A. Specifying the interaction 

Lacking any knowledge about the nature of dark energy and that of dark matter, there 
is no well-motivated choice for the interaction either. Our strategy here is as follows. 
Since we try to model a transient accelerating universe, we postulate a mathematically 
tractable interaction which presumably is compatible with such a dynamics. We solve this 
dynamics and clarify whether there exist parameter combinations that are favorable for a 
scenario of transient acceleration. It will turn out that for the two-parameter interactions 
considered below, there are indeed such ranges. Subsequently, we compare the result with 
the SNIa observational data. We shall consider couplings for which the function g(a) is 
represented by a combination of powers and exponentials of the scale factor. In particular, 
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we investigate the following two types of interaction terms: gi(a) = c\d n exp(— a/a) 
and #2(0) = C2d n exp(— a 2 /a 2 ), where n is a positive integer and a is a positive real 
number. Although these choices may seem arbitrary, they will admit a suitable parameter 
combination and they will provide us with a transparent picture concerning the potential 
role of interactions within the dark sector in a scenario of transient acceleration. 

In the first case, i.e., for gi(a) = c\a n exp(— a/a), assuming that 3w + n — 1 is a natural 
number, the energy density ( {TBI is given by 

,3„n— 3 



Kia^aJ 1 exp (— a/a) 

n+3w— 1 

3wK ia - 3(1+w) a 3 exp (-a/a) ^ 



i=0 



(n + 3w-l)\ 



a n+3w-i a i 



(19) 



where 



n+3w— 1 



J x 



p X0 +KiaQ exp (— a /a) + 3wKi exp(— a /a) 



i=0 



(n + 3w-l)' 



-a 



+3w—ii—3w 







(20) 



and K\ = c\p mo . From now on we shall put ao = 1. 
The expression f|T5|) specifies to 

^ + ^- + {l + 3w) Hx ° 



+ AnGw^exp^—a/a) 



a 3(l+w) 
n+3w— 1 



(1 + 3W) Y yU + 3V ! l)' a n + 3 W - V -3(l+ W ) +a n-3 
i=0 



(21) 



Let us now analyze the special case w — — 1 and n = 5 for which the sum in (1191) 
and ( 12"U|) has only two terms. Using p mo = p mo — K\ exp(— 1/a), the matter density Q 
reduces to 

p m = Pm a~ 3 + Kia' 3 [a 5 exp(-a/a) - exp(-l/cr)] , (22) 
and the density ( fT9l) of the x fluid becomes 



Px = P. 



> e JJ + 3K X exp (-a/a) (a 2 + aa - ^a 2 ^j 



with 



p e J f =p X0 -3K ie xp(-l/a) 



a A + a - 



(23) 



(24) 
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For a vanishing interaction the x component coincides with a cosmo logical constant. For 
the acceleration equation ( 12TI) one finds 



a 
a 



-Hi 



2Q e J f + 3^1 exp(-a/a) [a 2 - 2a 2 - 2aa] 



(25) 



wren j\i - ^jjiJXi ana ii xo - ^ p XQ . 

For an interaction #2(0) = c 2 a n exp(— a 2 /a 2 ) and w — — 1 and n = 5 it follows analo- 
gously that 

Pm = p mo a~ 3 + K 2 a~ 3 [a 5 exp(-a 2 /a 2 ) - exp(-l/(x 2 )] , (26) 

and 



P* = Pt! ~ K2 exp (-a 2 /a 2 ) [a 2 - \a 2 ^ 



where K 2 = c 2 p mo and 



PZ =P-o ~ exp(-l/(7 2 ) 



a 



(27) 



(28) 



Again, in the interaction- free limit K 2 — > we have p x — >■ p a . = const. The quantity p e J' 



can be seen as an effective cosmological constant which is re-normalized compared with 
the "bare" value, corresponding to p XQ , due to the presence of an interaction. The ratio 



- becomes 

a 

a - = -lm 



a 



n b0 -K 2 exp(-l/a 2 ) 



- 2Q e J o f + 3K 2 exp(-a 2 /a 2 ) [a 2 - a 2 ] \ (29) 



with K 2 = IZUK 9 and Q e JJ = ^Up e JJ. 



B. Decelerated and accelerated expansion 

To have a viable cosmological model, formulas ( )25l) and ( 1291) should admit a transition 
from - < to - > before the present time, i.e., for a < 1. If, moreover, the accelerated 
expansion is a transient phenomenon, there should be a change back from - > to - < 
at a future time, i.e., for a > 1. In both the expressions (125]) and (f29|) the a -3 terms on 
the right hand sides dominate for small values of a, i.e., there is decelerated expansion for 
a < 1 provided the conditions 

n mo +tob >K 1 exp(-l/<r) (30) 
8 



or 



Ono+ftfco >^ 2 exp(-l/ ( r 2 ) (31) 

are satisfied. These conditions put upper limits on the interaction strength. In the non- 
interacting limit they just express the positivity of the total matter energy density. 

Let's consider now the case a ^> 1. The dominating contributions in the braces on the 
right-hand sides of (j25]l and fl29|) are then given by the constant terms —2p e JJ. As long 
as p e Jf > we will have - > for a ^> 1, i.e., there is no transition back to decelerated 
expansion. This holds, in particular, in the non-interacting limit which reproduces the 
ACDM model. Then VL e JJ reduces to fl XQ , equivalent to Qa - F° r a ^ 1 this term will al- 
ways dominate the dynamics. An obvious way to obtain decelerated expansion for a 3> 1 
is to put p e JJ = both in (1251) and ( 1291) . This corresponds to a vanishing total cosmo- 
logical constant. In other words, part of the interaction cancels the "bare" cosmological 
constant, described by p xo . Under this condition it is exclusively the remaining part of 
the interaction which potentially can trigger a period of accelerated expansion. In such a 
case one obtains for the exponential interaction from ( Tl9l) 

1" 



Q xo = 3-fTi exp(-l/cr) 



<7 



a — 



while from (1271) for the Gaussian interaction 

Uzo = ^2 exp(-l/(J 2 ) 



—a 
2 



is valid. Adding up the energy contributions pb from 



(32) 
(33) 

p m from ( 1221 and p x with 



P. 



e JJ = from ffT9l). we find for the Hubble function 



H 2 1 - 3^1 exp(-l/a) (a 2 + a) 



H 2 



3K 1 exp(— a/a) (a 2 + aa) . 



(34) 



With p b from fl3]), p m from ( l2"B"j) and p x with p e JJ = from f[2"T|) . the corresponding 
expression for the second interaction is 



IP 
H 2 



1 - |a 2 ^ 2 exp(-l/ C T 2 ) 3 



o- 2 K 2 exp(— a 2 jo 2 ) . 



The acceleration equations ( 125]) and fl29|) simplify to 



a 
a 



1 ^2 f ^nio + ^feo _ Kiexp(-l/a) 
-rin 



+ 3iTi exp(-a/<r) [a 2 - 2a 2 - 2aa] 



(35) 



(36) 



and 



a 1 2 (Q mo +^b ~ K 2 exp(-l/a 2 ) 

— = —-tin 



„ .. + 3K 2 exp(-a 2 /a 2 ) [a 2 -a 2 ]) , (37) 

a 2 { cr J 

respectively. Moreover, with Q mo + f2fc = 1 — with Q xo from fl32l) or (133]) . their final 
forms are 



1 2 / 1 - 3K t exp(-l/a) [a 2 + a] 



and 



a 
a 



1 2 ri-|i^ 2( T 2 exp(-l/a 2 ) 
-rlr. 



3K X exp(-a/a) [a 2 - 2a 2 - 2aa] \ (38 



+ 3K 2 exp(-a 2 /a 2 ) [a 2 - a 2 ] 



respectively. To have decelerated expansion for a 1, 



K l exp(-l/a) [a 2 + a] < 



(39) 



(40) 



or 



K 2 a 2 exp(-l/(T 2 ) < - 



(41) 



has to be required. These conditions are equivalent to ( 13"U|) and (13"Ti) . respectively. The 
zeros of (1381) and (I39I) determine the values a q of a at which transitions between decelerated 
and accelerated expansion (or the reverse) occur, namely 

3i?i exp(-l/<r) [a 2 + a] + 3Ki a\ exp(-a q /<r) [2a 2 + 2ora 9 - a 2 ] = 1 (42) 



and 



\a 2 K 2 exp(-l/ ( x 2 ) + 3tf 2 a? exp(-a 2 / ( x 2 ) [a 2 - a 2 ] = 1 



(43) 



respectively. The conditions to have acceleration at the present epoch with a = 1 are 



iH 2 



> <^> Ki exp(— 1/er) 



1 



> 



1 



and 



aH 2 



> & K 2 exp(-l/a 2 ) 



a 2 - 



> 



(44) 



(45) 



If the inequalities ( l44"j) or fj45|) hold, we may have present acceleration under the condi- 
tion p e JJ = 0, i.e., a vanishing total cosmological constant. Obviously, the normalized 
interaction strengths K\ or K 2 have to be larger than a threshold value to realize this 
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configuration. The conditions (144j) or ()45l) are consistent with (1T81) if the latter is com- 
bined with either ( 1321) or ( 1331) . respectively. On the other hand, we have the upper limits 
fj4*0|) and f|4T|) . This means, there exists a range for admissible values of the interaction 
strength, determined by 

l/cr _ 1/cr 

<K,< - - ■ . (46) 



for the first case and 



9 (a 2 + a - |) ^ 3 (a 2 + a) 

2 e 1 ^ 2 -. 2e x / CT 



for the second case. 

It is instructive, to compare the relations ( 1581) and ( 1591) with the corresponding expres- 
sion 

HM^--4 <«> 

of the ACDM model. Obviously, the interaction terms in f|38|) and (159]) may play the role 
of f^A in ( I48j) . Now we know, that the ACDM model provides a fairly good description 
of the present universe, i.e., for a = 1. Therefore, we do not expect a radically different 
scenario around the present epoch, based on any alternative model. This suggest positive 
values of the interaction constants Ki and K 2 together with a > 1. It will turn out that 
these models are indeed favored. 



C. Direction of the energy transfer and structure of the source terms 

From it is obvious, that in an expanding universe (d > 0) the quantity -# should 
be positive to have an energy transfer from dark energy to dark matter. Positive values 
of Q are preferred on thermodynamical grounds j^. (But see H]-) For the case 
gi = C\a n exp(— a I a) we have 

da (a a) ® l ^ ^ 

It follows that for n = 5 and for c\ > 0, equivalent to K\ > 0, 

a < 5a & Q > . (50) 
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For g 2 = C2a n exp(— a jo ) 

dg 2 ( n 2a 
da \ a a 



--- 2 )92 (51) 



is valid. This means, for n = 5 and for c 2 > 0, equivalent to K 2 > 0, 

a 2 <%f 2 Q > • (52) 

Consequently, for sufficiently large values of the parameter a, positive interaction con- 
stants K\ and K 2 are both thermodynamically favored and provide a description of the 
present cosmological epoch which is similar to the well-tested ACDM model (see the 
comment following Eq. (1481) ). With w = — 1 our setting is reminiscent of models with 
a decaying cosmological "constant". Notice that for future values a > 1 for which the 
conditions ( 1501) or ( )52l) are violated, the interaction is no longer operative. 
The explicit forms of the interaction term Q in the balances ([2]) become 

Q x = K x a 2 H (5 - ^) exp (-a/a) (53) 

for the first interaction and 

Q 2 = K 2 a 2 H (h - exp {-a 2 /a 2 ) (54) 

for the second one. These terms may be written as 

Q 1 = A (a) H Px (55) 



with 



for the first interaction and as 



with 



Ma) = ^ 5 f—i (56) 

a? a 3 



Q 2 = 2 (a) Hp x (57) 



2 5_ 2— 

02 (a) = -,- 7 r-^ (58) 

6 ^ - 3 

respectively, for the second interaction. Interactions of the type Q = 0Hp x with an 
interaction constant have been frequently used to describe interactions between dark 
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matter and dark energy 20|, |28|, |29J . Here we have generalized this structure to a variable 
interaction parameter. For both the cases ( 1561) and ( 1581) this parameter vanishes for small 
values of the scale factor, i.e., the interaction is switched on during the cosmological 
evolution. In the late time limit both Qi and Q2 vanish because p x itself vanishes. 



D. Scalar field representation 

In this subsection we point out that there also exists an equivalent scalar field represen- 
tation for the dynamics of the interacting x component. Rewriting the balance equation 
for the x component (cf. Eq. (T5])) as 



p x + 3H(l +w eff )p x = 



(59) 



with 



W eff = W + 



Qi 



3H Px 



for the first interaction and introducing the correspondence 



(60) 



;i + W ef f) p x 



(61) 



we find 



and, via 



i 2 = I^a 2 (5 - ®) exp(-a/a) 



a 



v= ll-w eff 12 



an effective potential 



V = 3K ia 2 



21 + w eff 



a 2 ^ a 11 a 



a 2 a 18 18(7 



exp(— a/a) . 



(62) 
(63) 

(64) 



The kinetic part vanishes for very small and very large values of the scale factor. The 
potential starts with a constant initial value Vi ~ 3Ki<j 2 at a 1 and also tends to zero 
for a ^> 1. Eq. fl62|) is equivalent to 



% (5 - f ) exp (-a/a) 



1/2 



da 



H 
Ho 



(65) 
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with from (jM]) . The solution 0(a) of f )65|) together with f l64"|) yields an implicit repre- 
sentation of V (0). 

The corresponding relations for the second interaction are 



5-2^ exp(-a 2 / ff 2 ) 



and 



V 



-K,a 2 



a 2 2 a 2 
^2 + 9^ 
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exp(— a 2 /a 2 ) 



(66) 



(67) 



with a similar behavior as in the previous case. Eq. ( 166]) may be written 

1/2 

da 



^ 5-2^ exp(-a 2 /a 2 ) 



Ho 



(68) 



now with from fl35|) . 

Both potentials can be determined, as function of <f), numerically. The potentials can 
approximately be expressed by a finite series in terms of power and cosh functions. For 
example, the series can have the form 



V{4>) = a + 



ai< 



cosh( 



+ 



(>2 



cosh 



(69) 



For the first interaction we find ao = 0.121151, a\ = —0.527741 and a2 = 0.603898, while 
for the second interaction we have ao = 0.416204, ai = —0.655256 and 02 = —0.0962396. 
Figure (JT|) displays the numerical interpolation and the approximate analytical expression 
for each of the cases. 



IV. SUPERNOVA TYPE IA CONSTRAINTS 

To constraint the free parameters of the model, we use the gold sample of the SNIa 
data [23]. The crucial quantity to be evaluated is the luminosity distance 

r dz' 

Bi = c < 1 + *>/. 7FM- (70) 
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FIG. 1: Left panel: effective potential of the scalar field representation for the first interaction, right 
panel: effective potential of the scalar field representation for the second interaction. The continuous lines 
represent the numerical results, the dashed lines the interpolation in terms of power and cosh functions. 



where H is given either by fl34|) or by fl35|) . The quantity H = 100/i km/Mpc/s in these 
expressions is the Hubble parameter today. The observational data are expressed in terms 
of the moduli distance 

^o = logf^V 25 - ( 71 ) 



Mpc 

From this quantity the statistical function x 2 is constructed. It is defined by 

x2 = ^ (/4-/4) } (72) 

where /j^ is the predicted theoretical value for the moduli distance for the i th supernova, 
obtained by using the model developed above, /x^ is the corresponding observational data, 
(Ti being its observational error bar. With the help of \ 2 we construct the probability 
density function (PDF) 

P(x n ) = Ae- x2{Xn)/ \ (73) 

where x n denotes the set of free parameters of the model; A is a normalization constant. 
The general description of the Bayesian analysis can be found in reference 3p| . 

At first we obtain the results for the ACDM model which is given by the non-interacting 
limits K\ = and K2 = for the specific configurations considered here. The ACDM 
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model has only two free parameters: the density parameter tt x o = f^o an d the Hubble 
parameter h. Minimizing x 2 , we find Xmin = 1-128 with Q x q = ^ao = 0.651 and h = 
0.644. The two-dimensional PDF for these parameters is shown in figure [2] (left). After 
marginalization, we find the one- dimensional PDFs for Q\ Q and h, also displayed in figure 
[2] (center and right, respectively). After marginalization, the maximum probabilities occur 
for Q A0 = 0.651 and h = 0.619. 

In the interacting model there are generally four free parameters to be constrained: 
f2 X Q, K\ or K 2 , the Hubble parameter h and a. But our transient scenario implies the 
relations (|32|) or (|33l) between Q x0 and K\ or K 2 , respectively. This reduces the number 
of free parameter to three: K\ or K 2 , h and a. Integrating over one (two) of them we can 
obtain the two-dimensional (one-dimensional) PDFs for the remaining variable(s). From 
now on, we will restrict ourselves to the three-dimensional phase space case. 

For the first interaction (interaction strength parameter K\) the relevant equations 
are (134"]) and (138]) . The corresponding relations for the second interaction (interaction 
strength parameter K 2 ) are (135]) and (139]) . 

For the first interaction we obtain the two-dimensional and one-dimensional PDFs 
shown in figures [3] and HI The best fit is obtained for h = 0.643, a = 4.000 and K\ = 0.015 
with Xmin = 1-129, corresponding to Q x q = 0.683. The evolution of the acceleration 
parameter for this set of values is shown in the lower right panel of figure [3j It is obvious 
that after an intermediate accelerated expansion the Universe reenters a decelerating 
phase. After marginalization, the peaks of the one- dimensional PDFs occur for h = 0.670, 
a = 3.558 and Ki = 0.017. These values imply Q x0 = 0.599 and again a temporary period 
of accelerated expansion. The corresponding evolution of the acceleration parameter is 
displayed in the lower right panel of figure HI 

The results for the second interaction are very similar. In this case we find the two- 
dimensional and one-dimensional PDFs shown in figures [5] and [6j The best-fit values are 
h = 0.644, a = 5.000 and K 2 = 0.019, with xLn = 1-129, corresponding to Q x0 = 0.678. 
The evolution of the acceleration parameter is shown in the lower right panel of figure [5] 
After marginalization, the peaks of the one-dimensional PDFs are located at h = 0.670, 
a = 4.358 and K 2 = 0.022, implying Q x q = 0.563. The behavior of the acceleration 
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h ilAo h 

FIG. 2: The ACDM model. Left panel: the two-dimensional probability distribution for h and Oao- 
The brighter the color, the higher the probability. Center panel: one-dimensional probability distribution 
function for 51ao- Right panel: one-dimensional probability distribution function for h. 

parameter for this set of values is displayed in the lower right panel of figure El confirming 
again that the accelerated expansion of the universe is a transient phenomenon. We recall 
that all the cases of figures [3] - El correspond to an energy transfer from dark energy to 
dark matter. 



V. CONCLUSIONS 

In a homogeneous and isotropic cosmological background, a fluid with an equation 
of state w = — 1 is dynamically equivalent to a cosmological constant. According to the 
ACDM model, a small (in Planck units), but so far theoretically unexplained, cosmological 
constant is responsible for the observed accelerated expansion of our present Universe. 
We have demonstrated that an interaction between such type of dark energy and dark 
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FIG. 3: Two-dimensional probability distributions for different combinations of K\, a and h after 
marginalization over the remaining parameter. The lower right panel shows the evolution of the acceler- 
ation parameter for the best fit scenario in units of Hq . 

matter re-normalizes this constant. Even if the resulting effective cosmological constant 
is assumed to be zero, the time-dependent part of the interaction (the part that does 
not contribute to the re-normalization) was shown to be able to generate a phase of 
accelerated expansion of the Universe. We have studied two types of interactions for 
which cosmic acceleration is a transient phenomenon. While recent investigations favor 
a scenario in which the acceleration of the Universe is already slowing down today 
our toy model predicts the maximal acceleration to occur at a future time. The ACDM 
model is unable to provide any type of transient acceleration dynamics. A statistical 
analysis of the SNIa data shows that models of transient accelerated expansion may well 
be competitive with the standard cosmological constant model. While the mentioned 
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FIG. 4: One-dimensional probability distributions for K\, a and h after marginalization over the two 
remaining parameters. The lower right panel shows the evolution of the acceleration parameter for the 
highest peak probability values in units of . 



cancelation mechanism of the "bare" cosmological constant is purely phenomenological, 
we believe that it might possibly indicate a more fundamental feature. A vanishing total 
effective cosmological constant seems more appealing and easier to understand within an 
underlying basic theory than a small, non- vanishing value of this quantity. 

In general, our conclusion concerning the potentially transient character of the cosmic 
acceleration as the result of an interaction in the dark sector coincides with the numerical 



results of 22[. But we think that an analytic solution, even though it is the solution of 
a toy model, provides additional information about the specific role of the interaction in 
this scenario and thus gives rise to a physically more transparent picture. 
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FIG. 5: Two-dimensional probability distributions for the different combinations of K%, and h af- 
ter marginalization over the remaining parameter. The lower right panel shows the evolution of the 
acceleration parameter for the best fit scenario in units of Hq . 
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